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Abstract 



Weak first-order phase transitions proceed with percolation of new phase. The 
kinematics of this process is clarified from the point of view of subcritical bubbles. 
We examine the effect of small subcritical bubbles around a large domain of asym- 
metric phase by introducing an effective geometry. The percolation process can be 
understood as a perpetual growth of the large domain aided by the small subcritical 
^ . bubbles. 

Recently crucial effects of small rapid fluctuations -subcritical bubbles]]]] are actively 
^| investigated for the electroweak baryogenesis|| and for the inflationary cosmology. We 
could previously show, by studying the effect of subcritical bubbles, that the phase mixing 
is achieved in weakly first order phase transitions. However the global achievement of 
the phase transition -percolation- has not yet been explained from the point of view of 
subcritical bubbles. The difficulty lies on the previous treatment of subcritical bubbles; 
one always assumed spherically symmetry for the configuration of subcritical bubbles 
and interactions among them were neglected. However, it is possible to consider these 
effects easily within the point of view of subcritical bubbles. After the phase mixing is 
attained, there appear large and stable domains of asymmetric phase. If this domain were 
isolated in thermal fluctuations, then this domain would eventually collapse due to the 
surface tension of itself. However if we properly consider the effect of small fluctuations of 
subcritical bubbles around the large domain, this domain is stabilized against the collapse 
and eventually grows. This is the mechanism which we would like to clarify in this paper. 
This growth can be seen as the process of percolation of the system. Our study has been 
inspired by the interesting work by Gleiser et al. . However in our case, both the collapse 
and growth of the domain are automatically taken into account and we worry about the 
introducing an extra term which guarantees the decay as claimed in [H]. 

In this letter, we first show that a large asymmetric domain is stabilized by subcritical 
bubbles around the wall. For simplicity, we assume that the shape of the large domain 



is spherical with radius R(t) and volume V + (t). Hereafter we simply call the subcriti- 
cal bubbles as bubbles. The time variation of the volume V + (t) is composed from the 
kinematical and subcritical bubble contributions^. ; 

^ = + (T+AlOyW " (ToAOyWS, (2) 

where only the contribution from bubbles near wall of the large domain is included in the 
second and third terms in the right-hand-side of this equation. T + and T are the creation 
rates per unit volume of a bubble and an anti-bubble, respectively: 

T + ^m 4 + (T)ex V [-PF + ((R) + )] (3) 

and 

r ~m£(?)exp[-PF ((R) )], (4) 

where F + and (R)o,+ are the free energy and the averaged size of a bubble and an anti- 
bubble, respectively. As F ~ F + , m (T) ~ m + (T) and (R)+ — (R)o '■— (R) near T = T c 
at which two vacua degenerate, approximately T + = T =: T holds and therefore 

^«=4rff + r(AV-AV)f { R)*. (5) 

Furthermore, as AV^ — AV ~ 32ttR(R) 2 , this equation reduces to 

dV+(t) -dR 1287T 2 

This is the relation between the spherical symmetric volume and its radius modified by 
the bubbles around the large domain. A cute and elegant way to express this important 
relation is to introduce a fictitious geometry of non-Euclidean space whose metric is given 
by 

di 2 = dr 2 + a(t)r 2 dn 2 2 . (7) 

Deviation of the "scale factor" a from 1 represents the non-trivial effect from the last 
term in eq.(^). The time variation of the volume in this geometry (^) becomes 

dV+ d ( . /■«(*) , o , A d /4vr 



dt 



4nJ o drr 2 a(t)Uf t (^-R 3 (t)a(t) 



1 As we stated Gleiser et al. have considered that a domain without bubbles always collapses; in place 
of the first term on the right-hand side of this equation, they put a term —4irR 2 v where v is the wall 
velocity of the collapsing domain. We suspect this treatment; actual thermal fluctuations around a domain 
have both collapsing and expanding effects on the domain. Furthermore, they compare their results with 
computer simulations^). In this simulation the "final" fraction (7f) of the symmetric phase depends on 
the strength of the self-coupling A of the Higgs field which determines the mass. For sufficiently large 
A, this fraction becomes 0.5 within the Hubble time. However for smaller values of A, it never reach 0.5 
. In this latter case, the system is always in non-equilibrium state. On the other hand, they compare 
the thermal equilibrium value estimated through the Boltzmann equation with the above non- equilibrium 
value. The Boltzmann equation they use is given by 

d t F{cj>, R) ~ -vd R F{cj>, R) + (1 - 2/)G(0, R), (1) 

where F and G are the number density and the creation rate of the subcritical bubble inside the range 
((f) ~ <fi + 5(f>, R ~ R+5R), respectively. In the usual case, the first term of the right-hand side in the above 
equation vanishes and the equilibrium value should be 7 cq = 1/2. Of course the result in the simulation 
is <y slmul - = 1/2 if one continues the numerical calculation sufficiently beyond the Hubble time. 
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Equating eq.(H) and eq.(f|), one obtains an 'Einstein' equation for a 'scale' factor a(t); 



d /47T 



dt V 3 

The solution of this equation is 

327rr( J R) 5 



R 3 {t)[a(t) - 1] 



128tt' 



tr(rY 



a(t) = 1 + 



R 3 (t) 



f dt'R(t') = 1 + 7 x" 3 (r) f T dr'x(r') 
Jo Jo 



(9) 



(10) 



where 7 := 327rr(i?) 4 , x := R/(R) and r := t/(R). 

Next, we construct an equation of motion for 'matter' (R(t)). It is obvious that its 
Lagrangian is given by 

L = J d 3 x^C 

roc 

= Ana(t) I drr 2 C 
= a(t)L , 



where 



:n) 



'12) 



and M(R) := 1571 \J^ +R ' ■ Here we have used the form 0(x) = + exp[— |x| 2 /i? 2 (t)] for the 
large domain as previously^ f\ . Near the critical temperature the potential becomes 



V(R) = a{t)V {R) ~ -a(t)M(R). 

5 

Let us define a new variable:g(r) := J r dr'x{r'). Then the Lagrangian becomes 



C q. 



dq d 2 q s 
dr' dr 2 



M((R)Y 



1 + 7 



dq s 
dr 



-3 



1 dq ( d 2 q\ 2 
2~cfr [d^j 



2dq 
5 dr 



The corresponding Euler-Lagrange equation becomes 
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' dq \ d^q 
dr) d^ 



+ 107g 



d 2 q\ d 4 q 



dq\ 4 d 2 q d 3 q 



+ 20 , 
dr 2 I dr 4 \dr I dr 2 dr 3 



- 407*^2^+ 207 f^V^ 
dr dr 2 dr 3 \dr I dr 3 



d 2 q 

+ 247g^-|-127 



dr 



dr 2 



0. 



(13) 



(14) 



(15) 



2 For simplicity, we assumed that the configuration is Gaussian, but this is true if the domain has the 
thick wall. In general, the wall can be thin. However we expect that the non-Gaussian extension cannot 
give the drastic change on the present qualitative result. 
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It is hopeless to solve this equation analytically; numerical calculation is necessary. The 
results of numerical calculations are depicted in Fig. 1 and 2. They show the time 
evolution of x the radius R normalized by (R). In Fig. 1, we changed 7 fixing the initial 
velocity x(0). In Fig. 2, we fixed 7 and changed the initial velocity. ^From the figures 
one easily see that the domain always bounces by the fluctuation of bubbles. Actually 
we have always observed the bounce of the domain in a wide parameter range. Here, we 
note that we assumed d 2 x/dr 2 (0) = in these figures because we could not observe the 
drastic modification from the above qualitative results in another reasonable cases with 
d 2 x/dr 2 (0) 7^ 0. Therefore it seems that the bounce is a universal phenomena. Of course, 
the domain always collapses in the limit 7 — > 0. Our approximation is meaningful only 
for R > (R) when the domain is well distinguished from bubbles. Therefore we have 
restricted the parameter range which respects this constraint .0 

To understand the bounce of the domain, we concentrate on the behavior of the 
solution around the radius-minimum. We expand the 'scale factor' a(r) as a(r) ~ 1 + 

7rx -2 H , and neglect the velocity x(t) m 0. Then the Euler-Lagrange equation reduces 

to a simple form 

2 T 2x 1 

x + - =: x + d x V(x, t) ~ 0, (16) 

5 [x z + 7T x 

where V(x,t) = (2/5)ln[(x 2 +jt)/x\. Explicit time dependence in this equation prevents 
us to use energy conservation law for the analysis of kinematics of the domain. However, 
one can see rough behavior using the 'potential' V(x,t). The minimum of the potential 
is at x m (r) = (7t)^ 2 and the orbit of a domain is always bounded. The energy changing 
rate is given by 

™ ~ 2 -M((R))(32n(R) 3 T) { -^ ~ 2 -M((R))i^(AV - AV')T, (17) 

where (AV^ — AV^')r is the number of bubbles generated per unit time in the volume AV — 
AV. The factor (2/5)M((R)) is the free energy of one bubble. The factor ^{R) 3 /AV 
is the minimum fraction which can contribute to the energy of the domain. Thus we find 
that the above equation is quite reasonable! 

Though we have shown that a large domain always bounces and perpetually grows, 
it does not mean that the asymmetric domain eventually covers all the region of space. 
Actually the same argument is applicable for the domain composed of the symmetric 
phase instead of the asymmetric phase if the phase mixing is well established. The effect 
of small bubbles are the same and they stabilize both types of large domains. This domain 
growth can be regarded as the percolation process that a macroscopic order appears and 
grows in the system. Let us clarify this point further. 

We have considered the spherical symmetric domain in this letter. This analysis 
is also applicable for the evolution of the local curvature of a boundary between the 
symmetric and asymmetric phases instead of the radius a domain. The essential point 
in our argument was that the effect of bubble attachment is different in the convex side 
and in the concave side of the domain wall; the bubble attachment rate is higher in the 
convex side than that in the concave side. Especially for a small domain, bubbles attach 
almost only in the outside and this is the reason that all the domain bounces. Therefore 

3 We took the variable 7 as a parameter in this letter. In principle this is determined by the scale (R) 
of the bubbles as in 101 . 



4 



we can conclude that a large boundary wall always tend to become straight irrespect 
of the symmetry. This is the percolation process from the point of view of subcritical 
bubbles. We would like to report much quantitative results of percolation in our future 
publications based on the present argument. 
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Figure Captions 

Fig. 1. 

The variable x = Rj (R) vs time r = t/(R). We fixed the initial velocity x'{r) = — 1 
and changed the parameter 7: The parameter 7 is 0.016, 0.16, 1.6, and 16 in the order of 
increasing thickness of the curve. All domain solutions bounce and grows. 
Fig. 2. 

The same as Fig. 1, but we fixed 7 and changed the initial velocity x'(t) = — 1 for 
each figures. 7 = 0.16 for Fig. 2-a, 7 = 1.6 for Fig. 2-b, and 7 = 16 for Fig. 2-c. The 
parameter x'(t) is -1, -0.5, 0, 0.5, and 1 in the order of decreasing length of dashing. 
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